A quantum coherent approach to transport and noise in double barrier resonant 
diodes: shot noise a way to distinguish coherent from sequential tunneling 
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We implement a quantum approach which includes long range Coulomb interaction and investigate 
current voltage characteristics and shot noise in double barrier resonant diodes. Theory applies to 
the region of low applied voltages up to the region of the current peak and considers the wide 
temperature range from zero to room temperature. The shape of the current voltage characteristic 
is well reproduced and we confirm that even in the presence of Coulomb interaction shot noise can 
be suppressed with a Fano factor well below the value of 0.5. This feature is a signature of coherent 
tunnelling since the standard sequential tunnelling predicts in general a Fano factor equal to or 
greater than the value 0.5. This giant suppression is a consequence of Pauli principle as well as 
long range Coulomb interaction. The theory generalizes previous findings and is compared with 
experiments. 

PACS numbers: PACS numbers: 72.70.+m, 72.20.-i, 72.30+a, 73.23.Ad 
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I. INTRODUCTION 

Since its realizationi, the double barrier resonant diode 
(DBRD) proved to be an electron device of broad phys- 
ical interest because of its peculiar non Ohmic current 
voltage (I-V) characteristic. Indeed, after a strong super- 
Ohmic increase of current it exhibits a negative differen- 
tial conductance and eventually histeresis effects.^ Even 
the shot noise characteristics of DBRDs are of relevant 
interest in the sense that suppressed as well as enhanced 
shot noise with respect to the full Poissonian value has 
been observed (see Ref. for a review on the subject). 
These electrical and noise features are controlled by the 
mechanism of carrier tunneling through the double po- 
tential barriers. The microscopic interpretation of these 
features is found to admit a coherent or a sequential tun- 
neling approach'^. The coherent approach^ consists in 
considering the whole device as a single quantum system 
characterized by a transmission coefficient describing car- 
rier transport from one contact to the other. By contrast, 
the sequential approach- consists in considering tunnel- 
ing through the diode as a two step process where carri- 
ers first transit from one contact into the well, then from 
the well to the other contact. The intriguing feature of 
these two approaches is that from the existing literature 
it emerges that both of them are capable to explain the 
TV experimental data as well as most of the shot noise 



characteristics. Therefore, to our knowledge there is no 
way to distinguish between these two transport regimes 
and the natural question whether the tunnelling trans- 
port is coherent or sequential remains an unsolved one. 

The coherent approach to shot noise in DBRD has re- 
ceived wide attention since the first experimental evi- 
dence by Li et al^ of shot noise suppression with a min- 
imum value of the Fano factor 7 = Si/{2ql) = 0.5, here 
Si is the low frequency spectral density of current fluc- 
tuations and q the absolute value of the unit charge re- 
sponsible of current. Remarkably, most of the coherent 
approaches developed so fariSiSiiSiiiiiS predict a maxi- 
mum suppression 7 = 0.5 even if there is clear exper- 
imental evidence of suppression below this valuemd^iil 
down to values of 7 = 0.25M^ To this purpose, some 
authors obtained theoretical values of the Fano factor 
just below the value of 0.5, Ref. found 0.45 and 

Ref. 0.38, respectively. However, the physical inter- 
pretation of these results remains mostly qualitative and 
quoting Ref. this direction of research looks promis- 
ing but certainly requires more efforts. We remark that 
the theory of shot noise in DBRDs under the sequential 
approach^^d^d'''d^d^i^°i^^i^^i^^i^'^ provides a Fano factor 
of 0.5 as the minimum value of shot noise suppression. 

The aim of this paper is to develop a coherent theo- 
retical approach for current voltage and electronic noise 
in DBRDs accounting for Pauli principle and long range 
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Coulomb interaction going beyond existing models. To 
this purpose, we implem ent the quantum approach pro- 
posed in Ref . |22l25ll2fl | . We anticipate that the main re- 
sult of the present theoretical approach concerns with the 
prediction that suppression of shot noise with a Fano fac- 
tor below 0.5 is a signature of coherent tunneling against 
sequential tunneling. 

The content of the paper is organized as follows. Sec- 
tion 2 presents the theoretical model. Section 3 pro- 
vides the analytical expressions for the calculation of the 
current voltage characteristics in the low voltage region 
limited to the first peak of the current. Section 4 pro- 
vides the analytical expression for the electron noise cor- 
responding to the current voltage characteristics of Sect. 
3. Here, the Nyquist expression is recovered at vanish- 
ing applied voltage. At increasing voltages suppressed 
shot noise is found in the region preceding the current 
peak and enhanced shot noise at voltages just above the 
current peak. Section 5 reports a comparison between 
theory and existing experiments. Major conclusions are 
drawn in Sect. 6. 



II. MODEL 

The typical diode investigated here is the standard 
double well structure depicted in Fig. 1. We denote 
by r = (Fi + Tfi) the resonant states width, by the 
energy of the resonant level as measured from the center 
of the potential well and by F^^^; the partial widths due 
to the tunnelling through the left and the right barrier, 
respectively. We consider the case of coherent tunnelling 
in the presence of only one resonant state and we assume 
that the diode has contacts with unit square surfaces. 
For convenience, calculations are carried out using the 
cgs system. 

The kinetic model is developed by assuming that the 
electron distribution functions in the emitter and in 
the collector, /„, are equilibrium-like, but with different 
electro-chemical potentials Fa'- 

Us, F^) = \ . , (1) 



1 + exp 



here a = L stands for the left contact (the emitter), 
a — R for the right contact (the collector), fc^ is the 
Boltzmann constant, T the bath temperature and e the 
carrier energy. 

The full Hamiltonian of the structure is 



H — Hl + + Hres + Htun 



(2) 



here 



= Z^Z^(^(^^") +Pi/2™)ci(Pa,P±) Ca{Pa,P±) 

P± Pa 

are the Hamiltonians of the right and left contacts, 
and Ca the creation and annihilation operators of elec- 
trons in contact a, p± and Pa the momentum compo- 
nents perpendicular and parallel to the direction of the 



current, respectively, m is the effective electron mass in 
the conduction band, E{pa)+p'^/2m the electron energy 
in the a contact, 

Hres = ^{Er + p\ /2m)a^ [p i_) a{p±) 

p± 

is the Hamiltonian of the resonant states, Er = Sr — qu 
with u the voltage drop between the emitter and center 
of the quantum well, —q the electron charge, a"*" and a 
the creation and annihilation operators of electrons in 
the resonant level, and 

Htun^ ^ Taa+{pi_)y^Ca{pa,p±_) + hc\ 
P±,a \ Pa / 

the part of the Hamiltonian describing the electron tun- 
neling, Ta the amplitude of tunneling between the reso- 
nant state and the ath contact. 

Following Ref. [2^, the relation between F^ and 
is: Fq = 2'K\Ta\'^ Pa where pa is the one dimensional den- 
sity of states, and the electron operators a and c in the 
Heisenberg representation are given by: 

eicp{-i[E{pa) + p\/2m]t/h) 

a(t,P±) = Ta ^ , .^^^ Ca(Pa,p±) 



Ct,Pa 



E{pa)-Er+iT/2 



(3) 



Ca{t,Pa,P±) = Ca{Pa,P±) exp {-i[E{pa) +p\/2m]t/h) 



dTa{T,p±) exp [—i[E{p 

a 

)+pl/2m]{t-T)/h) 

_ J 

(4) 

The current operators for the left and right contacts 
Il,b^ and for the total current / are2^: 



iq 



{Taa+{t,p±)Cait,Pa,P±) - he) (5) 



/ = 77/L - (1 - 77)/^ 



(6) 



where 77 — u/V, V being the total appHed voltage (see 
Fig. 1). 

From Eqs. (3) to (5) we obtain the expressions for 
current operators similar to that obtained in Ref. 



P± 0,PH f'P-i 



^P'yiEp, E^)c'^ {pp ,pj_)c^{p^,p±) 



(7) 



where Ea = E{pa) and the full expression for 
A'^^{Ei3, Ej) is reported in the Appendix. From Eqs. 
(6) and (7) and the expression of Ap ^{Ep, E^) from the 
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Appendix, we find the usual expression for the average 
current 



< / >=< II > ~ 



< Ir >= 



qm 



/oc poo 
-QL JO 



(8) 

where e± — p^/2m is the kinetic energy of the transverse 
motion, e = (ez+ej_), gL (ul) the energy gap between the 
bottom of the conduction band and the first quantized 
level in the well before the left barrier (see Fig. 1), D{ez) 
the transparency of the double barrier given by: 



Die.) = 



In our model, we suppose that in the emitter there are 
no electron states with energy below —gL and that elec- 
tron states with energy higher than this value are three 
dimensional. 

To take into account Coulomb effects, we introduce the 
expression for the operator of the electron charge in the 
quantum well, that from Eq. (3) is found to be given by: 



QQw{t) = -q^a+{t,p±)a{t, 



P±) = 



'EEE 



T*T/3exp{i{Ea,~E[,)t/h) 
[Ea - Er- iT/2){Ep -Er+ iT/2) ' 



the emitter, the Poisson equation can be written in the 
form 



Anq 



/2 



(11) 



with 



NcFi,2 



Fl 



the electron concentration in the emitter, the effec- 
tive density of states of the conduction band, ^^1/2(2^) the 
Fermi-Dirac integral of index 1/2 n the static dielec- 
tric constant of the material. We note that the effect 
of size quantization on the electron distribution in the 
emitter is neglected. By integrating Eq. (11) and taking 
into account that ip (—00) = (—00) = 0, we find the 
relation between ip'j^ and ul on the left border of the left 
barrier 



NrF, 



3/2 



Fl + quL 



NrFr- 



3/2 



Fl 
knT 



knT 



nuL 



1/2 



(12) 

Here ^3/2 is the Fermi-Dirac integral of index 3/2iSl To 
simplify the task, we suppose that the barriers are un- 
doped and that the charge in the quantum well is placed 
at its middle plane, so that we can write 



+ v'l {dh + dqw /2) 



(13) 



and 



(9) 



v'l Q] [dR + dQw/2) (14) 



From Eq. (9) the average electron charge in the quantum 
well is found to be: 



< Qqw >— — 



qm 
2^ 



/OO POO 
dezD{ez) / deir-Vi(e,i^i 



qV 



dezD{ez 



de 



ri'fR{e,FR) 



(10) 



with dL, dfi, dqw the widths of, respectively, the left 
battier, the right barrier, and the quantum well; Q — 
[qNq^ — Qqw) the charge in the quantum well; Nq-^ 
the number of charged donors in the quantum well. Fur- 
thermore, we suppose that the electron concentration in 
the collector is the same of that in the emitter, hence 
Fr = {Fl — qV) and the Poisson equation in the collec- 
tor takes the form given by Eq. Ullfl with the change 
Fl Fr. By integrating this new equation with the 
condition ip (00) = V, f' {00) — 0, in analogy with Eq. 
(12) we obtain 



III. CURRENT VOLTAGE CHARACTERISTICS 



47r 



Vl 



k 

47r 



NcF 



3/2 



f Fl+ quR - qV 



knT 



The current voltage characteristic is determined from 
Eq. (8) once the dependence of D{ez) and 5l on the 
applied voltage is given. To calculate the transparency 
explicitly, we must find ul{V) and u(V) as functions of 
V . To this purpose, we consider the Poisson equation for 
the electrical potential Lp in the structure of Fig. 1. In 



'NcF^,2 



FL~qV \ 

keT / 



qn 
k^ 



{uR - V) 



1/2 



:0 (15) 



Equation (|15|l relates ur to V . We remark that it is 
more convenient to consider V and < / > as functions 
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of u because in this case they are single valued functions 
even for I-V characteristics of Z-type. Since the first and 
the third term in the left hand side of Eq. (|15|l are the 
charge of the emitter Ql and of the collector Qn, respec- 
tively, Eq. IjlSI) expresses the electroneutrality condition 
of the device. We note, that to derive Eq. (|15|l we have 
assumed that the width of the depletion region in the 
collector region (see Fig. 1) is smaller than that between 
the right barrier and the highly doped region in the col- 
lector. As typical in DBRDs, from both sides of the bar- 
riers there are spacers with doping level of the order of 
10"'^^ 10^^ cm~^ and with widths in the range between 
10 500 nmM^M.. If the low doped region in the collector 
is fully depleted, then Eq. must be substituted with 



V ^ur + ip'j^L 



2TrqnL'^ 



(16) 



where L is the width of the spacer in the collector, and 
(f'j^ = {(f'j^ — 'inQgw/K). In the derivation of Eq. (fTB|l . 
the voltage drop in the highly doped collector region is 
neglected. 



The charge fluctuation in the emitter and collector can 
be expressed through 6ui^ as: 



6Qe 



-CeSuL: 5Qc = -Cc6ui 



(21) 



where Ce,c = —dQe,cl duLji {V = const) are, respec- 
tively, the differential capacitance of the charge in the 
accumulation region of the emitter and in the deple- 
tion region of the collector whose expressions are de- 
tailed in the Appendix (we note that when Eq. (16) 
is valid Cc — k/{4ttL)). By taking in account that 
ip'j^ — —AnQe/ K, from Eq. (13) we find: 



SU^SUL{1 + Ce/Ci), Ci 



4T:{dL + dQw/2) 



(22) 



Analogously, by taking into account that ip'j^ — AttQ/k = 
AttQc/k, from Eq. (14) the relation between Su and Sur 
is found to be: 



6u^Sur{1 + Cc/C2), C2 = 



(23) 



4T:{dR + dQw/2) 
From Eqs. (21) to (23), Eq. (20) is rewritten in the form 



IV. NOISE 

To calculate the spectral density of current fluctua- 
tions at zero frequency under fixed voltage we use the 



expressioi 



,25. 



5/(0) = / dt< «(0)<5/(<) + 6I{t)SIiO) > (17) 

J —oo 

We anticipate that the total current fluctuation consists 
of two sources. The first comes from the fluctuation of 
population states in the contacts and the second from the 
fluctuation of the electron charge in the quantum wellS^. 
This last leads to fluctuations of the voltage drop between 
the emitter and the quantum well, Su. Accordingly, the 
operator of current fluctuations is given by>2^ 



Slit) = Shit) 



d<I > 
du 



Suit) 



(18) 



where Shit) is the current fluctuation operator under 
flxed u due to the population fluctuations in the contacts. 
Analogously we can introduce the operator of charge fluc- 
tuation 



SQqwit) = SQQW2it) - CqwSuit) 



(19) 



where Cqw is the differential capacitance of the quantum 
well 



Cqw — - 



d < Qqw > 
du 



From the condition of charge neutrality, the charge fluc- 
tuations in the emitter SQ^, in the collector SQ^, and in 
the quantum well SQqw satisfy: 



SQqw =SuiCL+CR), 



Ce,cCl: 



L,R 



iCe,c + C'1,2) 



(24) 



where Cl, iCn) is the capacitance between the emitter 
(collector) and the center of the quantum well, respec- 
tively. From Eq. (24) it is clear that Cl^r can be con- 
sidered as two capacitances in series Ci,2 and Ce,c- If 
we neglect the accumulation and depletion regions (i.e. 
Ce,c — * 00) and neglect the quantum well width (i.e. 
dow — > 0) then Cl.r coincide with those used in Ref. 
|2a |. From Eqs. (19) and (24), the relation between Su 
and SQqw2 is found to be: 



Suit) = 



SQQW2it) 

Cl + Cr + Cqw 



(25) 



From Eqs. (17), (18) and (25), 5/(0) is found to be given 
by the sum of three terms as: 

5/(0) = 5i + 52 + 53 

The first. Si, is the correlator of SI2; the second, 52, 
is proportional to the cross correlator between SI2 and 
SQqw2', the third, 53, is proportional to the correlator 
of SQqw2- They are given explicitely by: 



5i = 



52 = J 



dt < SI2 (0)«2 it) + SI2 it)Sl2 (0) > (26) 



dt < Sl2iO)SQQW2it) + Shit)SQQW2iO)- 



SQe + SQc + SQqw = (20) 



+SQQW2iO)5hit) + 5QQW2it)Sl2iO) > (27) 



5 



/OO 
dt < SQQW2{0)SQQW2{t)+ 
-OO 



-SQQW2{t)6QQW2{0) > 



(28) 



where 



(33) 

^3 = A^^/ dej de^D'x 



J = 



d<I> 



Cl + Cb. + Cqw du 



plays the role of a differential RC rate. 
By using the definitions for 6Ia2{t) 



-^fdl - h) + H%^[/l(1 - fn) + - /l)] J> + 



r2 



+A 



„2 /-OO <.oo pii 



(c^(P/3,P-L)c7(P7)Pi)- < c'l}ipf3,p±)Cj{p^,pi_) 



> 



for ^QQW2(t) 

SQQW2{t) = E E (ci(Pa,P±)c/3(p/3,?3±) 



- < ct{Pa,P±)cp{pi3,P±) >) X 



(29) 



T*T0 eMKEa - Ep)t/h) 



(30) 



{Ea - Er- iT/2){Ep - Er+ iT/2) 
and the propertjsSS 

< C^{Pa,P±)cp{pi3,P±)c:i;{p.y,p'j^)cs{ps,p'±) > - 
- < Ca{Pa,P±)Cf3{pp,p±) >< C+ {p^ , p'j_)cs{ps , p'±) > = 
= SasSp^psSi3j6pijp^dp^p'^fa{Ea + p']_/2m, Fa)x 

x{l-fpiEp+pl/2m,Fp)) (31) 
it is possible to find for S'i i = 1, 2, 3 the expressions 



Si = 



^ J de.J^ de^{D[fL{l-fR)- 



+ fB.{l-fL)]^D^{fL-fR)'} 



(32) 



Here we used the notation fL,R ~ Il^rIs, Fl,r) and 



A = 



hT 



1 



d<I > 



TlTr (Cl + Cr + Cqw) du 



where — oo<A<cxDisa dimensionless parameter de- 
scribing Coulomb interaction 

Equations (32)- (34) are the central result of this paper. 
We note that <S'/(0) does not depend on rj. When the 
A = 0, i.e. Coulomb interaction is negligible, Si{Q) — Si 
and the result of Ref. |23] is recovered. 

As an internal check we prove that 5/(0) satisfies 
Nyquist theorem^. Indeed, at zero applied voltage 
/l = Jr, from the expression (8) for the total current 
it follows that d < I > /du = 0, and therefore, A = 0. 
Accordingly, the differential conductance G at zero volt- 



age IS 
G 



We here used the property that for = dfR/dqV — 
fbil - fL)/{kBT). According to Eq. (32), Si at zero 
voltage is 



sm 



^4? r ^'^D r de^fL{l-fL) 

T^f^^ Jo Jo 



(36) 



Equations (35) and (36) imply 5/(0) = iksTG, which 
represents Nyquist theorem. 

Let us now show that 5/(0) — + 00 on the border of 
the instability region where d < I > /dV -^00. To this 
purpose, we note that dl/dV can be decomposed as: 



d<I>_d<I> d < I > du 
dV ~ dV ^ 

By writing the condition of charge neutrality as: 
—CeduL + dQqw + Cc{dV — duR) = 
by taking into account Eqs. (22), (23) and that 

dQqw = -Cgwdu + ^^'^ dV 



(37) 



(38) 



(39) 
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we obtain 



du 
dV 



= a 



dQqw \ 
dV J 



[Cl + Cr + Cqw) 



(40) 



By using Eq. (8) for < J >, we see that d < I > /Oul 
and d < I > /dV entering Eq. (37) are always finite. 
This impUes that d < I > /dV oo only if the sum 
(Cl + Cr + Cqw) — * 0. We remark that also the denom- 
inator of A is proportional to {Cl + Cr + Cqw)', thus, 
A, and in turn 5*7(0) go to infinity simultaneously with 
dl/dV. Note that Cl,r are always positive and Cqw be- 
comes negative when the resonant state approaches ql, 
that is when the number of electrons in the quantum well 
decreases at increasing u. 



r 



r2 



de{FL - e)D{e)^ 
(42) 

For the Fano factor, 7 = Si{0)/2q < I >, Eqs. (41) and 
(42) yield: 

7=1 j^^— < 1 - A A 

where 

A={f + 



r2 J B{f,o 

(43) 



l + (/ + 0' 



arctan(/ + £,) — 



V. RESULTS AND DISCUSSION 

In this section we present the theoretical results for 
the two relevant cases of zero and high temperatures (i.e. 
T ^ r which here corresponds to T > 77 K), the for- 
mer being appropriate to investigate the influence on shot 
noise of the Pauli principle and the latter of long range 
Coulomb interaction. Then, theoretical results are com- 
pared with experiments. 



A. Zero temperature. Three dimensional case 

We investigate the condition of high applied voltages, 
when qV > Fr, because more close to experiments. (We 
recall that typical magnitudes for the relevant parameters 
of DBRDs arc: for T less than a few meV, for Fl less 
than 100 meV, and for the voltage at the peak current 
around 0.5 V.) Because of the above, from Eq. (8) the 
expression for the current becomes: 



qm 
2^ 



de{FL - e)D{e) 



where 



(41) 



B{f, = 2(/ + e) [arctan(/ + ^ - arctan(0] 



In 



i + (/+e)^ 



and for convenience we define dimcnsionless chemical po- 
tential / and voltage drop ^ as 

^ _ 2{Fl + ql) ^^2{qu-er-gL) 



r ' ^ r 

The expression for the noise spectral density becomes 
•^^(0) = dsiFL - e)Die)[l - I?(e)] + 



■ arctan(^)] 



We note that A depends on / and through Cl,r, Cqw, 
d < I > /du which are given by Eq. (24) and by: 



Cqw — 



TT^h? r 



d <I> q^mVLTR 



du 



H{f,0 
(44) 



arctan(/ -I- ^) — arctan(^) 



dgL 



ATlTr dqu 



fD{~gL) 



If / ^ 1 and the resonant level is located below the the 
Fermi level of the emitter far from its borders Fl and 
i.e. C « -!,/ + ?» 1, then A{f,OlB{f,i) «l/2 
and Eq. (43) recovers the expression given in Ref. |2a | 



TI+TI + 2K{Tr-Tl) 



2A2 



with 



A 



Y2 



\TlVi. 



(45) 



As it will be shown below, when / ^ 1 the values taken 
by 7 in Eq. (45) are practically constant and correspond 
to the plateau exhibited by the dependence 7(/, C). On 
this plateau, 7 > 1/2 (the minimum is reached when 
A = (Fl - Fi^)/2). Note that on the plateau 



QW 



q^mVL 
7rft2r 



A; 



q^m 



1 



7rft2 Cl + Cr + C 



< 



QW 



Fl 



and the expression in the braces of Eq. (43) is positive. 
Therefore, on the plateau 1 > 7 > 1/2 and shot noise 
enhancement is impossible. 

Now we demonstrate that at voltage values for which 
the resonant level is close to the band edge of the emit- 
ter the Fano factor can drop below the value 1/2. To 
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this purpose, let us firstly consider the simplified model 
where we take k = 12.9 (GaAs), Cl = Cn — k/AttcI, 
d = b nm and neglect the term proportional to dgL/du 
in the expression for H{f,^). Figure 2 reports the de- 
pendencies of current and Fano factor on ^ for the two 
values / = 100 (Fig. 2 (a)), / = 10 (Fig. 2(b)) when 
Tl = Tfi and in the presence (continuous curves) or ab- 
sence (dotted curves) of Coulomb interaction. The fig- 
ure shows that for large value of the ratio between the 
Fermi energy and the resonant width (/ — 100), 7(^) 
exhibits a wide plateau region where 7 ~ 0.55 followed 
by a minimum with 7mm ~ 0.464. By further increasing 
the value of / (/ = 1000) the plateau region is found to 
widen and 7mm ~ 0.49. Note that, as it follows from Eq. 
(42), the Coulomb interaction always increases the noise 
if = Til With the decrease of / (see Fig. 2(b) where 
/ = 10), the plateau region becomes narrower and 7mm 
is found to decrease. The decrease of the value of 7miTi 
is due to two complementary reasons. The first is asso- 
ciated with the decreasing of the strength of Coulomb 
interaction. The second is associated with the increase 
of the effective barrier transparency near to the current 
peak and in turn with the further suppression of parti- 
tion noise. For the ideal case / ^ 1, at the peak current 
the transparency _D — > 1 and in turn 7min 0. 

Why does Coulomb interaction increase the shot noise 
if Fi = Tfi ? To answer this question, we analyze the 
noise contribution due to electrons with energies implying 
D — > 1. The part of corresponding to these electrons 
is zero, while the part of S3 is finite. We remind that 
the term 5*1 describes partition noise, thus it is propor- 
tional to D{1 — D) and for I? = 1 at zero frequency there 
is no noise associated with the fluctuation of 6l2- How- 
ever, even for D = 1 there is charge fluctuation in the 
quantum well due to the random character of electron 
escape from the well (random time-delay of charge) . The 
probability of this escape is proportional to F and this 
is the reason why 53 decreases with the increase of F. 
Therefore, Coulomb interaction enhances current noise 
due to the random time-delay of the charge in the quan- 
tum well, and which represents the quantum property of 
the electron motion in the RTD. 

The effects due to the asymmetry of the diode barriers 
when Tl = 0.25F and Tl = 0.75F are shown in Fig. 3 
which reports the dependencies of current and Fano fac- 
tor on ^ for the two values / = 100, / = 10. We note, 
that without Coulomb interaction (i.e. A = 0) these de- 
pendencies are the same for F^ — 0.25F and F^ = 0.75F. 
From this flgure it is clear that on the plateau Coulomb 
interaction decreases the noise for F^ = 0.25F while in- 
creases the noise for F^ = 0.75F. This asymmetry is due 
to the cross correlation term S2 which on the plateau is 
positive for Fl > Tjj and negative for F^ < Tr. The 
physical reason of this asymmetry stems from the follow- 
ing fact. If Fi > Tji, then the charge situated in the 
quantum well is characterized by an escaping probability 
to the emitter which is greater than that to the collector. 
By contrast, when F^ < Fjj the opposite happens. Of 



course, when F^ = Tn the escaping probability to the 
emitter and to the collector is the same and 52 = 0. We 
note the important role played by Pauli principle for the 
possibility of S2 to be positive. From Eq. (33) it is clear 
that if /l <C 1 and — then 6*2 < for any value 
of the ratio Fl/F/j. From Fig. 3 one can also see that 
7„„, = 0.457 < 0.5 for Fl = 0.75F and / = 100 (see Fig. 
3(a)) and for both values of F^ when / = 10 (see Fig. 
3(b)). 

We conclude, that for the simple model considered here 
both the increase of F and the decrease of Fl decreases 
7mm which value drops below 0.5 for F^ > F^. We 
further note, that there is bias when d < I > /du = 
and A changes of sign. Under this bias, for F^ — F/j 
the curves of the Fano factor with and without Coulomb 
interaction touch each other (see Fig. 2) while iorTi ^ 
Til they cross (see Fig. 3). 

To conflrm the possibility of evidencing the giant sup- 
pression 7mm < 0.5 in a real structure. Fig. 4 presents 
the calculations performed with a set of parameters re- 
lated to the experiments in Ref. 3 at T = 4.2 K and 
where both a symmetric (continuous curve) and an asym- 
metric (dashed curve) are considered. The structure pa- 
rameters are n = 2x 10^^ cm~^, dL = da — dgw — 5 nm, 
squared area of contacts 50 /xm^ and L = 50 nm, m — 
0.067 mo, with toq the free electron mass and k = 12.9. 
In the symmetric structure, the only two fltting param- 
eters are Fl = 0.5F = 0.48 meV and Sr = 104 meV. 
Their values control the location and the amplitude of 
the current peak, respectively, and are chosen by op- 
timizing the agreement between the experimental and 
calculated I-V characteristics at 77 K. For the asym- 
metric structure there are three parameters and we take 
Fl = 0.25F, F = 1.22 meV, Sr = 112 meV to preserve 
the location of the current peak at the same voltage. The 
function 5l(i*l) is calculated by solving the Schrodinger 
equation for a potential which: (i) for a; < follows from 
the Poisson equation without accounting for quantiza- 
tion effects; (ii) for a; > corresponds to the solid so- 
lution A^o. 42^00. 58^5 with zero electric field inside, as 
was done in Ref. j^l- For these values, the dependence 
<?l(ul) is found to be almost linear and well approxi- 
mated by: gL{uL) ~ 0.28{quL - llOfcsT) - 14^^^. The 
details of the solution of the Schrodinger equation are re- 
ported in the Appendix. Calculations give 7mm = 0.43 
in reasonable agreement with the value of 0.35 found in 
experiments* . From Fig. 4 we see that both the I-V and 
the noise characteristics are sensitive to the asymmetry 
of the structure, as expected. In any case, the main fea- 
tures of both characteristics are preserved. 

We conclude, that the main reason for shot noise sup- 
pression in RTDs at temperatures below about 4.2 K is 
essentially related to Pauli principle and, because of the 
coherent tunneling regime, near to the current peak the 
Fano factor can take values significantly lower than 0.5 
(giant shot-noise suppression). 
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B. Zero temperature. One dimensional case 

Here we shall discuss shot noise suppression for the case 
of one dimensional geometry. We recall that the possi- 
bility for 7mm < 0-5 when the resonance level touches 
the Fermi level in the absence of Coulomb interaction 
follows from the expression of the noise spectral density 
originally derived in Ref. |25l | and further discussed in 
Ref. 1^6] . At zero temperature we can write the expres- 
sions for the current and current noise spectral density 
in the following forms (spin degeneracy included) 



SiiO) 



deD{e) 



deD{e) 



(46) 



1 ^ ^L-'^R yi ^L^R 



deD{ef 



(47) 



By using Eqs. (46) and (47) and after integration, for 
7we obtain the expression: 



here 



Y2 



j Si(/,a 

(48) 



+ arctan(/ + $) — arctan(^)] 



Bi{f; = arctan(/ -(- ^ — arctan(^) 



The dependencies Cqw and d < I > /du on /, ^ are 
Cqw 



7rr2 



with 



1 



du 



1 



If we neglect the Coulomb interaction (i.e. A = 0), intro- 
duce 2; = / -I- ^ and take ^ ^ — oo (the so called infinite 
band model'^^), then Eq. (47) recovers the expression for 
the noise spectral density derived in Ref. • We note, 
that similarly to the three dimensional case there is a 
plateau region in the behavior 7(^) when / ^ 1. How- 
ever, in the present one dimensional case d < I > /du is 
small and hence Coulomb interaction weakly influences 
the current noise. In contrast to the three dimensional 
case, here the Coulomb interaction influences the noise 



when the resonant level is: either close to the Fermi level 
in the emitter Fj^, or to the conduction band edge — 
We note, that here the expressions for the capacitances 
derived in the three dimensional case are no longer valid, 
and furthermore we have assumed Cl,r. to be constant. 
Since when / 3> 1 the minimum value of Cqw is approx- 
imatively given by —Aq'^TL/irT'^, for the appearance of a 
current instability due to Coulomb interaction it should 
be (Ci + Cfl) >4g2ri/7rr2. 

Figure 5 reports the calculated dependencies of the cur- 
rent and Fano factor versus ^ for the symmetrical struc- 
ture Tl = Ffi when / = 100 (Fig. 5(a)) and / = 10 (Fig. 
5(b)), respectively. Here, we take 7rr(CL — Cii)/2q'^ — 
0.9, which corresponds to the absence of current insta- 
bility, and neglect terms proportional to dgh/dqu in the 
expression for Hi{f,^). The I-V characteristics exhibits 
a symmetric shape with respect to^ = — //2asa conse- 
quence of a r which is independent of the applied volt- 
age and of the property of the one-dimensional density 
of states. From this figure it is evident that jmin < 0.5 
even when / = 100. We note, that due to the reasons 
discussed above, the Coulomb interaction increases the 
noise for the symmetric case when F^ = F^j. In addition, 
the reason of shot noise enhancement is the Coulomb in- 
teraction. Accordingly, the enhanced result obtained in 
Ref. p6| in the absence of Coulomb interaction is not 
physically plausible. From Fig. 5 it is clear, that the 
Coulomb interaction in the one dimensional case influ- 
ences the noise when the resonant level is close either to 
the Fl (e - -/) or to -gL (? ~ 0). 

The effects due to the asymmetries of the structure are 
reported in Fig. 6 which shows the dependencies of the 
current and Fano factor versus ^ when F^ = 0.25r and 
Tl = 0.75r with the same (Cl -I- Cr) value as in Fig. 5. 
We note, that in the absence of Coulomb interaction, the 
Fano factor is the same for F^ = 0.25r and Tl ~ 0.75r. 
From Fig. 6 one can see that the Coulomb interaction 
decreases the value of "/min- In particular, when F^ < 
r_R, 7mm is located at ^ -/, while, when Fl > F^, 
at ^ ~ 0. We note, that for £, ^ —f the cross correlator 
< SI2SQQW2 > is negative when F^ < F^; and positive 
when Fi > Ffl. In both cases A is positive. On the 
contrary, for ,^ ^ A is negative because d < I > /du < 
0. 

We conclude, that also for a one dimensional RTD the 
Fano factor can drop below the value 0.5 in the presence 
or in the absence of Coulomb interaction. 



C. High temperatures 

We now discuss the case of high temperatures when 
ksT ^ r, which in the present case refers to T > 77 K. 
Again, we consider applied voltages high enough to ne- 
glect the contribution to the current and noise of the elec- 
tron fiux moving from the collector to the emitter. Since 
the energy scale for the change of il'(£)is significantly less 
than that of /^(e), all the integrals in the expressions for 
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the current and the noise spectral density can be calcu- 
lated analytically. To perform the calculations, first of 
all we note that 



1 + exp 



Fl - £2 



and A < T /Tl in the stable region where (Cl + Cr + 
Cqw) > 0. ^ 

We now estimate the possible minimum value for 7. 
From Eq. (55) it is clear that to obtain 7mi„ it is neces- 
fcBr$L(£2)sary that 



(49) 



A 



r 

2rl 



1 



2rz 
r 



1 



<&(/?) 



(57) 



1 + exp 



knT 



We note, that this value of A lies in the range of possible 
values characterizing the stable region. By substituting 
Eq. (57) into Eq. (55) we find 



and 



d£^D{e^)<^L.2L{£z 



L.2L 



= fcsT$2L(e. 
(50) 

(51) 



ah 



de,Die,)HL,2Liez) 



9L 



r3 



L,2L 



here /? is the maximum value between (e 



(52) 

-qu) and -ql, 



|-arctan(0, D2{0 ^ - 



1 



By using Eqs. from (49) to (52), we can write the fol- 
lowing expressions, respectively for the current, the noise 
spectral density and the Fano factor 



SiiO) = 



2q^mkBTTLTR 



(53) 



r2 



2T, 



^2l{P)+\[^l{P)--^^2l{P) 



^*2l(/3) 



(54) 



2VlTrD2{0 



-A' 



2^L 



$(/3) + A 1 



2T, 



.<i>(/3) 



(55) 



here $(/3) = $2l(/3)/'&l(/?)- Analogously, we derive the 
expressions for Cqw and d < I > /du. Since in this case 
d < I > /du ~ TrCqw /ti, we obtain 



A 



rc, 



QW 



(56) 



In 



D2{0 

2T DiiO 



(58) 



Now we analyze two limiting cases. 

The first is the case when the occupation factor of the 
state with energy /3 is much smaller than unity. Then 
$(/3) « and: 



= 1 



Tr D2{0 

2r DiiC) 



(59) 



Since the maximum value of D2{£,) / F>i{^) equals 1.217 
(^ = ^0 = -0.802), from Eq. (59) we find that 7„i„ < 0.5 
when Tr > 0.82ir, and that the expression (57) holds 
near C = Co- 

The second limiting case is when the probability of 
occupation of an electron state with energy (3 is close to 
unity and <&(/?) « 1. In this case 



7n 



= 1 - 



^2(0 

2D^{0 



(60) 



and 7mm < 0.5 when the expression (57) holds in any 
point where < 0, since there Z?2(C)/^i(C) > 1- 

In concluding this subsection, we emphasize that, for 
RTD theory predicts values of the Fano factor below 
0.5 also at high temperatures. We note that, provided 
$(/3) = /L(/3)//L(/3)^ Eq. (55) is valid also for the one 
dimensional case and thus again "/min < 0.5 is confirmed 
to be possible. 



D. Comparison of theory with experiments 

We compare theory with two sets of experiments per- 
formed on DBRDs with barriers sufficiently narrow to 
expect that coherent tunneling is of importance. The 
first set refers to pioneer experiments of Brown^ which 
are detailed at 77 K with indications at 300 K. The sec- 
ond set refers to recent experiments at 300 K reported 
in Ref. _32 33] . In both cases the comparison is limited 
to the voltage region up to the peak current since theory 
neglects energy levels in the quantum well higher than 
the first one. 

Figure 7 reports the comparison between experiments 
performed in Ref. |^ and present calculations at 77 K. 



10 



Numerical results make use of the same parameters in 
Fig. 4 for the symmetric structure. For the used values, 
the dependence is found to be almost linear and 

well approximated by: ~ 0.44(gui — LSfc^T) — 

0.07fcsT for 77 K. 

From Fig. 7 (a) we see that present calculations well re- 
produce the I-V characteristic including the current peak. 
From Fig. 7 (b) we see that the calculated Fano factor re- 
produces both the suppression and enhanced behaviors. 
However, its minimum value is of 0.65 against the exper- 
imental value of 0.35. By choosing appropriate values of 
Tl^h the theoretical model can be forced to fit the experi- 
mental Fano factor at the expense of overestimating the I- 
V characteristics for about one order of magnitude. This 
result corrects previous findings of the same authors^, 
where the contribution of charge fluctuations to the total 
noise was underestimated with respect to the present ap- 
proach. In this context, we note that Ref. [l3| presented 
a theoretical calculations of the same experimental at 
77 K. The results of these calculations were found to be 
in excellent agreement with experimental data except for 
the region of instability. On the borders of this region 
the noise tends to infinity and, as a consequence, there 
were no measurements of noise in this region. However, 
contrary to such an experimental evidence, the theoreti- 
cal calculations.^- provided finite values of the noise and 
the the absence of the instability region, which makes the 
theoretical approach at least suspect. 

To provide a physical insight of the physical reason 
for shot-noise suppression. Fig. 8 reports the results of 
the calculations associated with the presence (continu- 
ous curves) and the absence (dotted lines) of Coulomb 
interaction. From Fig. 8 (a) we see that at 77 K the 
Coulomb part contributes to suppress shot noise in the 
whole region of applied voltage, and that the Pauli part 
becomes noticeable near the current peak. Furthermore, 
while Pauli contribution leads systematically to suppres- 
sion, the Coulomb contribution becomes responsible of 
enhanced shot noise at voltages near to 0.57 V where the 
instability region is approached. From Fig. 8 (b) we see 
that at 4.2 K the predominance of Pauli over Coulomb 
interaction in suppressing shot noise is confirmed. Again, 
Coulomb effects are found to be responsible of enhanced 
shot noise around above the current peak, in agreement 
with experiments^. 

To complete the analysis of this structure at 77 K, we 
have investigated the role played by the electron concen- 
tration and the width of the quantum well. Accordingly, 
Fig. 9 reports the current voltage characteristic and the 
Fano factor for the same structure of Fig. 7 but with 
two different values of the electron concentration n in 
the injector. For the first structure we use ghiuh) ~ 
0.37(guL - ^ksT) - l.l^ksT (n = 2 x 10^^ cm.-^) and 
for the second gUuL) ~ 0.22(guL - 4fcsT) - O.OgfcsT 
(n = 2 X 10^^ CTO^^). From the figure it is clear that 
the increase of n leads to: (i) the increase of the value of 
the current peak, (ii) the widening of the region of shot 
noise suppression and, (iii) the decrease of the voltage 



value corresponding to the current peak and the mini- 
mum of the Fano factor. We note that for values of n 
which differ over two orders of magnitude the peak cur- 
rent changes only for a factor of about 2.5. It means 
that, due to accumulation, the electron concentration on 
the border of the left barrier differs only for a factor of 
2.5 for voltages corresponding to the current peak. Fig- 
ure 10 reports the current voltage characteristic and the 
Fano factor for the same structure of Fig. 7 but with 
two different values of the quantum well width dqw- 
We have taken into account that when dgw decreases 
both r and increases by using for the first structure 
Er — 145 meV, Tl — O.bT = 1 meV and for the second 
structure Er = 45 meV, Tl = O.bT = 0.05 meV. This 
leads to the growth of the current peak and its shift to 
higher voltages. From the same figure we see that simul- 
taneously with the growth of the current peak also the 
width of the suppression region rises and shifts to higher 
voltages. Furthermore, the minimum value of the Fano 
factor is found to become more pronounced at increas- 
ing carrier concentration and decreasing the width of the 
quantum well. 

Figure 11 reports the current voltage characteristic and 
the Fano factor for the same structure of Fig. 7 at 300 K. 
Here the solution of the Schrodinger equation provides: 
gL{uL) = 0.53{quL - l.5kBT) + Q.lksT. Calculations 
show, that even by increasing the temperature, the main 
features of transport and noise already shown at 77 K 
are preserved. However, the current peak and the min- 
imum of the Fano factor become less pronounced at in- 
creasing the temperature. These trends are in agreement 
with experimental results^ which claim a reduction of 
the peak-current value of about 1 mA and an increase of 
the minimum of the Fano factor when going from 77 to 
300 K . However, even at 300 K present calculations give 
a Fano Factor for about a factor of two greater than that 
found by experiments. 

A recent set of experiments-^'^3: carried out at 300 
K on a DBRD with barriers thinner than those of Ref. 

, thus more adequate to check coherent tunnelling at 
high temperature, is reported in Fig. 12 together with 
theoretical calculations. The structure consisted of two 
2 nm AlAs layers separated by 6 nm InGaAs quantum 
.^gjj32j36 Measurements were carried out at 300 K using 
a Noise Figure Meter (XK5-49), that allows to measure 
simultaneously noise figure and power gain of two-port 
networks in the 50 f2 feed circuit. Simultaneously with 
the noise the I-V curve was measured. The diode was 
mounted in the break of a microstrip line, with one elec- 
trode been grounded, and another one bonded to the ends 
of a microstrip line. Noise measurements at frequencies 
60 and 200 MHz showed the same results within an ex- 
perimental uncertainty at worst of 20%, thus indicating 
that 1// noise contribution is negligible. Numerical re- 
sults makes use of the following values for the parameters 
entering the model: = 87 meV , r^; = Fl = 1.08 meV , 
n = 5 X 10^^ cm~^ , = (Ir — 2 nm, dgw — 6 nm 
L = 50 nm, gL{uL) = 0.445(guL - LbkeT) + GAblkeT, 
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m = 0.045mo. Also here the only two fitting parame- 
ters are and F — O-Sr^^^, all other parameters being 
provided by the experimental conditions, Figure 12(a) re- 
ports the I-V characteristic which shows a region of pos- 
itive differential conductance (pdc) up to about 0.7 V 
followed by an instability region. In the same pdc re- 
gion, the Fano factor is found to exhibit a suppression 
with a minimum value of about 0.4 at around 0.65 V 
[see Fig. 12(b)]. As for the case at 77 K, from Fig. 
12 (a) we see that present calculations well reproduce 
the I-V characteristic including the current peak. From 
Fig. 12 (b) we see that the calculated Fano factor re- 
produces both the suppression and enhanced behaviors. 
However, the minimum value of the Fano factor is found 
to be of 0.75 against the experimental value of 0.40. By 
choosing larger values of T^^ji the theoretical model can 
be forced to fit the Fano factor at the expenses of over- 
estimating the I-V characteristic for about one order of 
magnitude. This result corrects previous findings of the 
same authors^^, where the contribution of charge fluctu- 
ations to the total noise was underestimated. Thus, the 
comparison between theory and experiments at temper- 
atures above about 77 K provides qualitative agreement 
but is not able to explain the drop of the Fano factor 
below 0.5 found in experiments. 



VI. CONCLUSIONS 

We have implemented a quantum mechanical approach 
to investigate DBRDs transport and noise characteristics 
within the coherent tunneling regime that includes both 
Pauli principle and long range Coulomb interaction. The 
expression for the current voltage and noise character- 
istics generalize previous finding£^^''S^*2&. In agreement 
with expectations, at increasing voltages theory predicts 
a current characteristic which exhibits a peak followed 
by an instability region and that before the current peak 
shot noise is suppressed because of Pauli principle and/or 
Coulomb interaction. In addition, theory confirms shot 
noise enhancement well above the full Poissonian value 
at the current peak as a consequence of the positive feed- 
back between tunneling and space charge. 

At zero temperature, the suppression of shot noise 
starts in concomitance with the sharp increase of the 
current associated with the alignment within the value 
of r of the Fermi level in the emitter with the resonant 
level in the quantum well. Accordingly, the Fano factor is 
found to exhibit a minimum near the current peak. Re- 
markably, the value of this minimum can be significantly 
below the value 0.5 of the full Poissonian value. This gi- 
ant suppression of shot noise can be taken as a signature 
of coherent tunneling since sequential tunneling can pre- 
dict suppression but at most with a Fano factor not less 
than 0.5. At 4.2 K, for a realistic DBRD structure we 
find a minumum value of the Fano factor of 0.44 that is 
in agreement with experiments^. At temperatures above 
about 77 K , we have found that coherent tunneling still 



predict that shot noise can be suppressed well below the 
value of 0.5. This giant suppression has been evidenced 
by experiments performed at 77 and 300 K. However, 
the present theory is not able to explain this noise fea- 
ture with the same set of parameter values able to explain 
the I-V characteristics. Probably the model used is still 
too simple, and further efforts are needed to provide a 
better interpretation of experiments. 

In any case, to our opinion, the main result obtained 
here is the fact that in RTDs the drop of the Fano fac- 
tor below 0.5 can be taken as a signature of coherent 
against sequential tunneling. Therefore, below we briefiy 
discuss the physical reason why shot noise suppression is 
expected to be more effective for coherent than for se- 
quential tunneling. Starting from the fact that the two 
mechanisms responsible for shot noise suppression are 
Pauli principle and Coulomb interaction, we note that 
both act simultaneously for coherent and sequential tun- 
neling. Let us consider the first mechanism, which is the 
most relevant at zero temperature, in the case when the 
Fermi energy is sufficiently small so that all the electrons 
exhibit the same transparency. Then, coherent trans- 
port admits a transparency near equal to unity, which 
implies 7 = (1 — D) ~ according to Lesovik findings^. 
For sequential transport the total transparency is always 
less than unity because of the finite value of the differ- 
ential rates controlling the relaxation of carrier number 
fluctuations inside the two terminal device through the 
contacts. As a consequence, under coherent transport 
for the possible case of full transparency (i.e. D = 1) 
there is no noise. By contrast, under sequential trans- 
port the presence of scattering inside the quantum well 
always introduces noise. This example illustrates why 
Pauli principle is more efficient in suppressing shot noise 
under coherent than sequential transport. By consider- 
ing Coulomb interaction, which is more relevant at high 
temperatures, we recall that in the absence of collisions 
it provides giant shot noise suppression^! as in the vac- 
uum tubes'^^ because electron refiection in this case is due 
only to Coulomb interaction. It is clear that the presence 
of scattering provides additional random mechanisms for 
electrons returning to the emitter and, therefore, pro- 
vides an additional source of noise. Even this example 
shows that Coulomb interaction is more efficient in sup- 
pressing shot noise under coherent than sequential trans- 
port. We finally want to stress that the main reason of 
the difference between these approaches stems from the 
fact that the sequential tunneling is based on a master 
equation- -'^'^^ for treating fluctuations of carrier numbers 
inside the quantum well, while coherent tunneling uses 
the quantum partition noise as only source of interaction 
and fluctuations. The master equation describes implic- 
itly a sequential mechanism for a carrier entering/exiting 
from the well and, as a consequence, its intrinsic limit co- 
incides with that of two independent resistors (or vacuum 
diodes) connected in series and each of them exhibiting 
full shot noise. This system yields a maximum suppres- 
sion of shot noise down to the value of 0.5. By contrast. 
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partition noise, inherent to a quantum coherent formal- 
ism, can be fully suppressed down to zero in the presence 
of a fully transparent barrier and weak Coulomb interac- 
tion like in vacuum diodes. 
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VII. APPENDIX 

Here wc detail the calculations to evaluate the current 
operator, the differential capacitances of the DBRDs and 
the solution of theSchrddinger equation for the voltage 
dependence of the first electron level in the emitter. 



A. Current operators 

By using Eqs. (3), (4), (5) and (7), the explicit expres- 
sion for A'^^^{Efj, Ej) is the following 



ir/2 



' E^ - Er + iT/2 



{Al) 



where 



Ef3 - Er + iT/2 
From the above, we note the following properties: 
A£^(i?, E)pL = A«^(i;, E)pR = -A'^^iE, E)pR 
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and 



= -Af^{E,E)pL = D{E) 



Aij,{E,E)A^i^{E,E)pLPR 



(A3) 



A^j,{E, E)A^^{E, E)pLpn = D{E){1 - D{E)) (A4) 

Eqs. (A3) and (A4) are useful for current and noise cal- 
culations. 



B. Capacitances 

By recalling that Ql = — «</?z,/47r from Eq. (12) we 
have 



Fl + quL 



— n 



Fl + quL 



NcF^/2 



knT knT 



nUL 



-1/2 



[Bl) 

If the electron accumulation is relevant (i.e. qul > ksT) 
then we can write 



Ce{UL) ~ q 



«iVci^3/2 (^) 



{B2) 



By substituting for Qn the value given in Eq. (15), for 
Cr we find 



Ce = q^n-N,Fy2 ( 



Fl + quL - qV 
ksT 



SirksT 



here y = ^' is the electron wave fimction and e the 
electron energy. Since d(p/dx in the emitter is function 
of (fi [see Eqs. (11) and (12)], instead of x it is convenient 
to use (f asvariable quantity so that 



dy_d(p _^ 2 
d(f dx 



2m 



{qip + s) = 



(C2) 



where 



dx 



SnkBT 



NrF, 



3/2 



Fl + qf 
keT 



NrF 



3/2 



\kBT 



ksT 



1/2 



(C3) 



The potential in the emitter takes values in the range 
from far from the barriers [x — oo) to ul on the bor- 
der of the left barrier. We note that far from the barriers 
the potential falls exponentially tp ~ cxp(a;/Ai3) when 
qip <§; fc^T, where A/j is the Dcbyc screening length. For 
localized state the electron wave function ~ exp(A;ia;) 
for x — > — oo, where k\ = \/—2me/h and, thus, we 
have as boundary condition y{p = 0) = fci. Since we 
suppose that the field in the barrieris absent, the elec- 
tron wave function in the barrier is exp(— ^22;) 
where k2 = \j2m (AEc — qul — e)/h and AEc is the 
conduction band offset on the barrier border. By using 
the conditions of continuity of the wave function and of 
its derivative, we obtain the second boundary condition 
y{ip = Ul) = —k2- Equations (C2) and (C3) together 
with the boundary conditions allow us to find the en- 
ergy of the first electron level in the potential well of the 
emitter by numerical calculations. 



A^ci^3/2 



Fl + quR - qV 



NcFs/2 



FL-qV 
keT 



qn 
"k^ 



{UR - V) 



-1/2 



If q(y — ur)/ (kBT) > 1, then we have the usual expres- 
sion for the capacitance of the depletion region 



qnn 



Sn{V - ur) 



(B4) 



C. Energy of the first electron level in the emitter 

The one dimensional Schrodinger equation for an elec- 
tron moving in the potential —q(p{x) can be written in 
the following form 

dy 9 2m , s „ 
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FIG. 1: Sketch of the band profile of the double barrier struc- 
ture considered here The bottom of the conduction band in 
the emitter in the well and in the collector coincides at V — 
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FIG. 2: Dependence of the Fano factor and of the cur- 
rent on the electrical potential in dimensionless units ^ — 
2{qu -£r- gL)/r with / = 100 (a) and / = 10 (b) for the 
symmetrical structure Tl = Tr. Here / = 2{Fl + ql)/^ and 
lo = qraV bJ {2'jt'^ h^) . Continuous (dashed) curves corre- 
spond to the presence (absence) of Coulomb interaction. 
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FIG. 3; Dependence of the Fano factor and of the current 
on the electrical potential for asymmetric structures with / = 
100 (a) and / = 10 (b). The dimensionless units are the same 
of Fig. 2. 
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FIG. 4: Dependence of the current and Fano factor on the 
appUed voltage for the structure of Ref. [A] at T = 4.2 K. 
Values of Fl.a and Sr are chosen from the fitting of the current 
voltage characteristic at 77 K. 
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FIG. 5: Dependence of the Famo factor and current on the 
electrical potential for a one-dimensional symmetric struc- 
ture Fl = Tr with / = 100 (a), and / = 10 (b). Contin- 
uos (dashed) curves correspond to the presence (absence) of 
Coulomb interaction. Here 7i = 2qTLTii/{-whT). 
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FIG. 6: Dependence of the Famo factor and current on the 
electrical potential for asymmetrical one dimensional struc- 
tures with / = 100 (a), and / = 10 (b) 
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FIG. 7: Calculated (solid) and measured (dashed) dependen- 
cies of current and Fano factor on the applied voltage for the 
structure o^ at 77 K. The parameters are the same used for 
Fig. 4 in the case of the symmetric structure. 




FIG. 8: Calculated dependencies of the Fano factor with (con- 
tinuous curves) and without (dashed curves) Coulomb inter- 
action at 4.2 and 77 K. Other parameters are as in Fig. 7. 
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(b) 
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FIG. 9: Dependencies of current and Fano factor on applied 
voltage for structures with n = 2 x 10^^ cm~^ (solid) and 
n = 2 X 10^^ cm~^ (dashed). Other parameters are as in Fig. 
7. 
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FIG. 10: Dependencies of current and Fano factor on applied 
voltage for structures with dgw = 3 nm (dashed curve) and 
dQw = 8 nm (solid curve). For the first structure wc used 
Sr = 145 meV, Fl = O.SF = 1 meV. For the second struc- 
ture we used Sr = 45 meV, Fl = 0.5F = 0.05 meV. Other 
parameters are as in Fig. 7. 
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FIG. 11: Dependencies of current and Famo factor on applied 
voltage at T = 300 K. Other parameters are as in Fig. 7. 
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FIG. 12: Experimental and calculated dependencies of cur 
rent (a) and Fano f actor (b) on applied voltage for the struc 
ture in Ref. |32l3,'j | at T = 300 if. 



